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The dark knight falters 
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ABSTRACT 

Tentative line emission at 111 and 12 9 GeV from 16 unassocia ted Fermi-h AT point 
sources has been reported recently by ISu fc Finkbeinerl (|2012ct ) . Together with sim- 
ilar features seen by Fermi in a region near the Galactic Centre, the evidence has 
been interpreted as the spectral signature of dark matter annihilation or internal 
bremsstrahlung. Through a combination of supervised machine-learning algorithms 
and archival multiwavelength observations we find that 14 out of the 16 unassociated 
sources showing the line emission in the Su & Finkbeiner sample are most likely ac- 
tive galactic nuclei (AGN). Based on this new evidence, one must widen the range of 
possible solutions for the 100-140 GeV excess to include a very distinct astrophysical 
explanation. While we cannot rule out a dark matter origin for the line emission in 
the Galactic Centre, we posit that if the detection in the Su & Finkbeiner sample is 
indeed real it might be related to accretion, bubble, or jet activity in nearby [z < 0.2) 
AGN. Alternatively, given the right conditions, the similarity could be due to a chance 
occurrence caused by extragalactic background light (EBL) absorption. Or else one 
must concede that the features are an artefact of instrumental or calibration issues. 
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1 INTRODUCTION 

Frantic activity has ensued over the past few months fol- 
lowing the report of an excess of Fermi gamma-ray events 
clustered around 100 and 140 GeV in a region near th e 
Galactic Centre (|Wenigerl 120121 : ISu fc Finkbeinerl l2012bl ). 
Dark matter annihilation and internal bremsstrahlung 
have rapidly emerged as poss ible explanations 
l|Bringmann e't all [ioTa IWenigerl 120121) Alternative inter- 
preta tions have also been advanced JProfump fc Linden I 
20121; iBovarskv. Malvshev fc Ruchavskiv I 120121 : 



Aharonian. Khangulvan fc Malvshevl 20121) More 



cently things have heated up even further with a tantalising 
claim of similar line emission at 111 and 129 GeV in 16 
unassociated sources detected in the Second Fermi-LAT 
catalogue (2FGL). The detection could provide independent 
support for a dark matter origin for the line emission see n 
near the Galactic Centre region (ISu fc Finkbeinerl l2012d ). 
Certainly, such coincidence might not only help us unlock 
the mysteries of dark matter, but it w ould also prove the 
existence of dark matter subhaloes (|Klypin et al.l 1 19991 : 
iMoore et allll999l ; ISpringel et al.ll2008l ). 

In the absence of obvious flaws in the analysis, the col- 
lected evidence has risen as a sort of dark knight - albeit an 
indirect one - that might finally grant us non-gravitational 
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access to dark matter. Intrigued by this possibility, we ex- 
plo re the nature of th e 16 Fe rmi unassociated sources listed 
by ISu fc Finkbeinerl (|2012<J >. Based on machine-learning 
classification algorithms and multiwavelength examination, 
we show that 14 out of the 16 unassociated Fermi sources 
displaying the lines are likely gamma-ray AGN. Therefore, 
rather than strengthening the argument, the detection of an 
identical signal in the Su & Finkbeiner sample appears to 
disprove a dark matter origin for the Fermi features unless 
a set of very unique astrophysical conditions are met. 

The paper is structured as follows. In Section [2] we 
explain the machine-learning classifier Sibyl. In Section [3] 
we compile clas s prediction for the 16 u nassociated Fermi 
sources listed in ISu fc Finkbeinerl i|2012d) . Section 0] details 
multiwavelength searches for the potential counterparts of 
these 16 objects. Finally, we provide some interpretation in 
Section [5] 



2 SIBYL 

Initially, we want to examine whether the set of unasso- 
ciated sources showing line emission at 111 and 129 GeV 
|Su fc Finkbeinerll2012cl ) is distinct in any way when com- 
pared to the bulk of associated Fermi sources. Confirmation 
of a truly unique type of gamma-ray source would hint that 
we may have finally found the much sought-after dark mat- 
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ter subhaloes predicted by numerical simulations. As our 
base for such comparison, we use Sibyl, a Random Forests 
classifier that generates predictions on class memberships 
for unassociated Fermi-LAT sources u sing gamma-ray spec - 
tral features extracted from the 2FGL (|Mirabal et al, l2012h . 
Random Forests (RF) create an ensemble of classifiers with 
a tree structure, where the splits captures the com plexity of 
the fe ature space among the set of training objects l|Breimanl 
2001). To tag a new object, RF lets each classifier vote and 
then outputs a prediction based on the majority of votes 
(P> 0.5). RF also computes proximities between pairs of 
objects and quantifies which variables are instrumental to 
individual classification. Previously we performed a simi- 
lar analysis for una ssociated 2FGL sources at |fo| ^ 10° 
l|Mirabal et al.ll2012T ). Here, we extend the coverage down to 
\b\ 5° in order to encompass the entire Su & Finkbeiner 
sample. The classifier present ed has been implemen ted with 
the R ran domForest package (Liaw fc W iener 20021). 

As in iMirabal et al.l (|2012T ). we trained Sibyl using 800 
labelled AGNs (BL Lacs and flat-spectrum radio quasars 
only) and 108 p ulsars from the the complete Fermi LAT 
2FGL catalogue (|Abdo et al.|[20lb1 ; lAckermann et al.|[201ll ; 
lAbdo et al.ll2012r ). There are additional gamma-ray classes 
in the 2FGL, b u t sinc e the 16 unassociated sources in 
ISu fc Finkbeinerl (|2012ch lie at \b\ ^ 5° we do not expect 
noticeable contamination from Galactic plane sources. The 
main tas k thus is to find o u t whet her the unassociated sam- 
ple from lSu fc Finkbeinerl (|2012cT ) falls into these two cate- 
gories or it is clearly different from these types of objects. 

During training and testing with the 908 labelled 
sources, we used a total of 7 spectral features: Index, 
Curve, Vari ability, and Flux R atios (FR12, FR23, FR34, 
and FR45) l|Mirabal et al.l [20121 ). Assuming class bimodal- 
ity, Sibyl achieves an accuracy rate of 97.1% based on ma- 
jority voting (97.7% for AGNs and 96.5% for pulsars). In- 
spection of the results shows that misidentifications tend 
to occur when less than 70% of the individual classifiers 
(P< 0.7) agree on a particular prediction. Therefore we set 
this threshold as our internal limit for a valid prediction. 




log(Flux1000) 

Figure 1. Kernel density plot of the 1—100 GeV photon flux 
FluxlOOO (photons cm -2 s -1 ) for the 573 unassociated sources 
listed in the 2FGL (solid). Comparison with the 16 unassociated 
sources in the Su & Finkbeiner sample (dashed). There is no 
obvious difference between the peaks of the samples. 



Furthermore, none o f the objects appears to be an out- 
lier (|Mirabal et al.l [2012) or has been strongly suspected 
to be a potential dark matter subhalo candidate in any of 
the p r evious searches conducted to date jBucklev fc Hooperl 
20101; iNieto et all [20lH: iBelikov. Hooper fc Bucklevl 12011 



Zechlin et al.l 2011 ; Ackermann et al.l 20121 ; Mirabal et al l 
20121 ). Lastly, the photon flux distribution for the 16 unas- 
sociated sources is compared to the overall distribution for 
the full Fermi-LAT unassociated sample (573 sources) in 
Figure [T] Except for a slight mismatch at the very bright 
end, we find no obvious selection effects that could produce 
line emission in this particular set. 



3 APPLICATION TO THE SU & FINKBEINER 
SAMPLE 

For each of the 16 una ssociated Fermi sources listed in 
ISu fc Finkbeiner (2012c), Sibyl provides a prediction that 
the object is an AGN (Pagn) or a pulsar (Pp u i 3a r) based 
on individual votes tallied from the classifiers. We adopt a 
threshold P > 0.7 to accept a prediction i.e., at least 70% of 
the trees agree on the final decision. Conservatively, sources 
failing to reach the threshold remain formally unassociated 
and most likely constitute interesting gamma-ray sources. 
In total, Sibyl predicts 14 objects in the Su & Finkbeiner 
sample to be AGN. The resulting predictions and percent- 
ages of voting agreements are listed in Table [T] Only two 
objects 2FGL J1716.6-0526c and 2FGL J1721.5-0718c re- 
main without a firm prediction. We note that both sources 
are also the only ones fitted with LogParabola functions in 
the Su & Finkbeiner sample. But such a pair is not uncom- 
mon as there are at least 163 associated Fermi sources with 
LogParabola best fittings in the 2FGL i ncluding numerou s 
AGN, pulsars, and supernova remnants (|Abdo et alj|20l3 h 



4 MULTIWAVELENGTH EXAMINATION 

To the untrained eye, it might seem overly simplistic to rely 
on machine-learning algorithms to make a class prediction 
for a particular source. We simply refer the reader to the vast 
amount of research and applications of classification meth- 
ods that have managed to reach tr emendous accuracies in a 
variety of astrophys ical subfields (|Bloom fc Richards! 1201 ll ; 
iRichards et al.ll2012h . However, one must never forget that a 
smart computer generated guess is no substitute for observa- 
tion Q. We take this recommendation at heart and move the 
association process even further by searching for the actual 
counterparts in archival multiwavelength observations that 
have partially or fully covered the Fermi 95% confidence er- 
ror ellipses of the 16 unassociated sources. Fo r this purpose, 
we employ a set of well- validated strategies ([Mirabal et alj 
|2000| ; iReimer et al1l200ll ; iMukheriee et al.ll2002l ).~ 

The wide distribution of sources over the sky results 
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Table 1. Predictions and voting percentages for the Su & Finkbeiner sample, ordered by RA. 



Source 


Pagn 


f Pulsar 


Prediction 


2FGL 


J0341.8+3148c 


0.786 


0.214 


AGN 


2FGL 


J0526.6+2248 


0.926 


0.074 


AGN 


2FGL 


J0555.9-4348 


0.992 


0.008 


AGN 


2FGL 


J0600.9+3839 


0.952 


0.048 


AGN 


2FGL 


J1240.6-7151 


1.000 


0.000 


AGN 


2FGL 


J1324.4-5411 


0.974 


0.026 


AGN 


2FGL 


J1335.3-4058 


0.864 


0.136 


AGN 


2FGL 


J1601. 1-4220 


0.996 


0.004 


AGN 


2FGL 


J1639.7-5504 


0.942 


0.058 


AGN 


2FGL 


J1716.6-0526c 


0.612 


0.388 




2FGL 


J1721.5-0718c 


0.498 


0.502 




2FGL 


J1730.8+5427 


0.986 


0.014 


AGN 


2FGL 


J1844.3+1548 


0.996 


0.004 


AGN 


2FGL 


J2004.6+7004 


0.996 


0.004 


AGN 


2FGL 


J2115.4+1213 


0.996 


0.004 


AGN 


2FGL 


J2351.6-7558 


0.960 


0.040 


AGN 



on a hodgepodge of radio and X-ray observations from var- 
ious existing catalogues. For radio counterparts, we relied 
on measur ements from the Gre en Bank (GB6) catalogue at 
4.85 GHz ijGregorv et al.ll 19961 ). t he 4.85 GHZ Parkes-MIT - 
NRAO (PMN) survey catalo gue (iGriffith fc Wrightl 1 19951 ) . 
the 1.4 GHz NVSS catalogue (ICondon et al.ll 19981) . and the 
843 MHz SUMSS catalogue i|Mauch et al.ll2003l ). In total, 
we find that 13 out of the the 16 unassociated sources have 
prominent potential radio counterparts within their Fermi 
95% confidence error ellipses. 



We complement the radio matches with observations 
from the most ambitious X-ray pr ogram for counter - 
part identification currently underway ( Falc one et al.ll20Tll ) , 
which aims to image the totality of unassociated Fermi 
sources with the SwiffX-recy telescope. To date, nine sources 
have been imaged by Swift with times ranging from 1.1 to 
19.1 ks of useful exposure. Source extraction to identify all 
significant X-ray sources within the Fermi error ellipses was 
performed with wavdetect. Source positions and positional 
errors were derived using xrtcentroid. X-ray counts (0.1-2.4 
keV) were extracted from a circular region with a 20 pixel 
radius (47"). The background was extracted from an annu- 
lus with a 20 pixel (inner radius) to 30 pixel (outer radius) 
around the source. Throughout, we used XSELECT to in- 
clude counts with grades 0-12. Six out of the nine Swift fields 
have single potential counterparts within their Fermi 95% 
confidence error ellipses. 

The ROSAT All-Sky Survey Faint Source Catalogue 
i|Voges et all l2000h adds four more potential single X-ray 
counterparts to the final count. Table [2] summarises the 
counterpart candidates in each case. Of the 16 sources, seven 
have both radio and X-ray tentative counterparts. Figure 
[2] shows X-ray flux versus radio flux density of associated 
Fermi AGN. Superposed are the potential seven with simul- 
taneous radio and X-ray counterparts. The results are in 
line with radio and X-ray counterpart flux levels expected 
for typical Fermi AGN. But we must emphasise that without 
dedicated spectral classification in the optical these must be 
considered solid but tentative counterparts at this stage. 
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Figure 2. X-ray flux versus versus 1.4 GHz flux density 
{Si.AGHz)- Small dots represent associated Fermi AGN. The 
black squares mark the seven unassociated sources from the Su &: 
Finkbeiner sample with tentative counterparts in both radio and 
X-rays. 

5 INTERPRETATION AND CONCLUSIONS 

We have presented class predictions of the Random Forest 
classifier Sibyl for 16 unassociated Fermi sources showing 
line emission at 111 and 129 GeV. We find that 14 out of 
16 unassociated sources in the Su & Finkbeiner sample are 
AGN candidates with prediction accuracy rates greater than 
96.5%. In addition, we have detected 10 X-ray and 13 radio 
potential counterparts distributed over the 16 unassociated 
Fermi 95% confidence error ellipses that would be consistent 
with the AGN predictions. We emphasise the word potential 
here as a more exhaustive detective work must be completed 
to confirm the appropriate counterpart for each unassociated 
source. 

It was postulated that the gamma-ray lines among 
the unassociated were perhaps connecte d to dark matter 
subhaloes dragged into the Galactic disc |Su fc Finkbeiner! 
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Table 2. Potential radio and X-ray counterparts for the Su & Finkbeiner sample. 



Source X-ray Counts s" 1 (0.1-2.4 keV) Radio S„(mjy) 



2FGL 


J0341.8+3148c 














NVSS J034213±314739 


Si. 4GHz 


= 23 


2FGL 


J0526.6+2248 














NVSS J052643±225337 


SlAGHz - 


= 107 


2FGL 


J0555.9-4348 














PMN J0555-4345 


S4.85GHz 


= 61 


2FGL 


J0600.9+3839 


Swift J0601.04 


-3838 


(2.7 ± 1.0) 


X 


io- 


3 


GB6 J0601±3838 


S4.85GHz 


= 322 


2FGL 


J1240.6-7125 


Swift J1240.4- 


7148 


(2.6 ±0.1) 


X 


io- 


1 


MGPS J124021-714901 


Sg43MHz 


= 18 


2FGL 


J1324. 4-5411 


1RXS J132455.7- 


-542020 


(5.4 ± 2.0) 


X 


io- 


2 


PMN J1325-5419 


S4.85G.Hz 


= 56 


2FGL 


J1335.3-4058 














SUMSS J133603-405758 
SUMSS J133544-410113 
SUMSS J133535-405407 


Ss43MHz 
S&43MHz 
Sg43MHz 


= 154 
= 21 
= 18 


2FGL 


J1601. 1-4220 














PMN J1600-4227 
PMN J1600-4217 


S4.85GHz 
S4.85GHz 


= 48 
=46 


2FGL 


J1639.7-5504 


1RXS J164023.6- 


-550259 


(2.0 ±0.9) 


X 


io- 


2 








2FGL 


J1716.6-0526c 


1RXS J171657.0- 


-053418 


(3.0 ± 1.2) 


X 


io- 


2 








2FGL 


J1721.5-0718c 


1RXS J172147.4- 


-071923 


(2.3 ±0.9) 


X 


io- 


2 








2FGL 


J1730.8+5427 














GB6 J1731±5429 


S4.85G.Hz 


= 19 


2FGL 


J1844.3+1548 


Swift J1844.44 


-1546 


(7.7 ±0.5) 


X 


io- 


2 


GB6 J1844+1547 


S4.85GHz 


= 76 


2FGL 


J2004.6+7004 


Swift J2005.1- 


-7004 


(7.9 ±0.5) 


X 


io- 


2 


NVSS J200506±700440 


SlAGHz 


=7 


2FGL 


J2115.4+1213 


Swift J2115.4+1218 


(1.9 ±0.5) 


X 


io- 


2 


NVSS J211522+121802 


SlAGHz 


=16 


2FGL 


J2351. 6-7558 


Swift J2351.3- 


7600 


(1.3 ±0.3) 


X 


io- 


2 


PMN J2351-7559 


S4.85G.Hz 


=47 



2012c). However, assuming an isotropic distribution, at least 
160 Fermi AGN are expected at |fe| ^ 10°. To date only 
abou t 100 are accounted for in the 2FGL l|Ackermann et all 
l201ll ). Thus, it makes astrophysical sense that AGN are 
making up an important fraction of the Su & Finkbeiner 
sample even at relatively low Galactic latitudes. 

In light of these results, the dark matter ori- 
gin for the narrow gamma-ray features observed by 
Fermi is in question. Were the s e dark matter sub- 
haloes ([Baltz. Taylor fc Wai I 120071 ; iDiemand et al. I 120081 ; 
iKuhlen. Madau fc Si 13 l2009h . coincidence between the 
Galactic Centre and the Su & Finkbeiner sample 
woul d certainly confirm a dark matter particle ori- 
gin (|Hooper fc Lindenl I2012T ). However, the interpretation 
changes dramatically if the unassociated sources showing 
an identical line signature are AGN, as implied by both 
machine-learning classifiers and the multiwavelength argu- 
ments just presented. Dark matter could be fed into AGN 
jets and the Galactic Centre, but such an explanation feels 
contrived given the hadr onic and leptonic dominan ce in the 
gamma-ray photon field l|Hinton fc Hofmannll2009l ). 

Instead, a distinct astrophysical mechanism unrelated 
to dark matter annihilation and linked to nearby AGN 
(z < .2 to avoid redshifted lines) such as a ccretion, 
bubble dSu et al.l l2010l: iProfumo fc Linden I [20121 ) . or jet 
(|Su fc Finkbeinei 2012af ) phenomenology would appear to 



be more logical. However, we note that although many Fermi 
AGN display photons above ~ 10 GeV, only a handful of soft 
AGN (r > 2) exhibit maxi mum photon energies greater 
than 100 GeV at z > 0.5 jAckermann et al.ll201ll ). Con- 
sequently, ISu fc Finkbeinei] (|2012rJ ) might be detecting a 
fiendish cluster of events imprinted by EBL absorption in 
the same energy band, but completely unrelated in origin to 
the emission observed near the Galactic C entre region. 

Od dly enough, the lines reported by ISu fc Finkbeiner! 
|2012d ) appear to be only present collectively in unasso- 
ciated sources and do not appear as pronounced among 
associated sources, inclu ding well-known gamma-ray AGN 
jSu fc Finkbeinerll2012d ). Therefore, we must also admit the 



possibility that the spectral signatures detected by Fermi 
origin ate from confounding ins trumental or calibration prob- 
lems (|Hooper fc Lindenll2012D . The Fermi calibration team 
will have the final word on the matter very soon, but in- 
dependent efforts must be made to scan the public Fermi 
archive for gamma-ray lines among individual AGN at z < 
0.2, as well as in diffuse emission outside the Galactic plane. 

We shall hear more about this energy region by 
the end of the yea r with t he recently unve i led H .E.S.S. 
II (jBecherini et alj |2012| ; iBergstrom et all |2012| ). and 
even more sensitive observations will be available later 
on after completion o f the Cherenkov Telescope Array 
jCTA Consortiumll201lh . In the future, a dark knight might 
rise again. Until then, we eagerly await for the final chapter 
of this intriguing saga. 
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